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SUMMARY 


The  Karman-Pohlhausen  method  is  extended  primarily  to  sixth-degree 
velocity  profiles  for  determining  the  characteristics  of  the  compressible 
laminar  boundary  layer  over  an  adiabatic  wall  in  the  presence  of  an  axial 
pressure  gradient.  It  is  assumed  that  the  Prandtl  number  is  unity  and 
that  the  coefficient  of  viscosity  varies  linearly  with  the  temperature. 

A  general  approximate  solution  which  permits  a  rapid  determination  of 
the  boundary -layer  characteristics  for  any  given  free-stream  Mach  number 
and  given  velocity  distribution  at  the  outer  edge  of  the  boundary  layer 
is  obtained.  Numerical  examples  Indicate  that  this  solution  will  in  ,, 
practice  lead  to  results  of  satisfactory  accuracy,  including  the  critical 
Reynolds  number  for  stability.  For  the  special  purpose  of  calculating 
the  location  of  the  separation  point  in  an  adverse  pressure  gradient,  a 
short  and  simple  method,  based  on  the  use  of  a  seventh-degree  velocity 
profile,  is  derived.  The  numerical  example  given  here  indicates  that 
this  method  should  in  practice  lead  to  s\iff iciently  accurate  results. 

For  the  special  case  of  flow  near  a  fooTward  stagnation  point  it  is  shown 
that  the  Karman-Pohlhausen  method  with  the  usual  fourth-degree  profiles 
leads  to  results  of  adequate  accuracy,  even  for  the  critical  Reynolds 
number. 


INTRODUCTION 


In  reference  1  it  was  concluded  that  from  the  viewpoint  of  both 
accuracy  and  convenience  of  calculation  a  suitable  method  for  deter¬ 
mining  the  characteristics  of  a  compressible  laminar  boundary  layer  is 
that  based  on  an  extension  of  the  Karman-Pohlhausen  integral  method  to 
velocity  profiles  of  higher  degree  than  the  fourth,  especially  sixth 
degree.  An  ordinary  differential  equation  for  general  types  of  flow  was 
derived,  but  only  the  flow  over  a  flat  plate  at  zero  incidence  was 
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investigated  in  detail.  The  purpose  of  the  present  investigation  is  to 
apply  explicitly  this  method  to  flows  with  axial  pressure  gradients. 

The  ordinary  differential  equation  derived  in  reference  1  is  con¬ 
verted  here  into  a  convenient  nondimens ional  form,  and  a  general  approxi¬ 
mate  solution  of  this  equation  in  closed  form  is  then  derived.  By  means 
of  this  solution  the  physically  significant  houndary-layer  characteristics 
of  the  flow  over  an  insulated  wall  can  he  calculated  fairly  easily  and 
quickly  for  any  given  free-stream  Mach  number  and  for  any  given  velocity 
distribution  outside  of  the  boundary  layer.  For  particular  accuracy  in 
determination  of  the  separation  point  in  an  adverse  pressure  gradient, 
a  simple  method  based  on  the  use  of  a  seventh-degree  velocity  profile, 
which,  in  accordance  with  a  suggestion  of  Timman  (reference  2),  is  made 
to  satisfy  an  additional  condition  at  the  separation  point,  is  derived. 

For  illustrative  purposes,  two  simple  but  basically  different  types 
of  flow  are  treated  in  detail.  In  the  first  example  flow  with  a  linearly 
diminishing  velocity  at  the  outer  edge  of  the  boundary  layer  is  con¬ 
sidered.  Such  a  flow  is  of  particular  interest  here  because  it  repre¬ 
sents  the  simplest  case  of  an  adverse  pressure  gradient  and  because  the 
usual  application  of  the  Karman-Pohlhausen  method  with  fourth-degree 
velocity  profiles  has  been  known  to  lead  to  highly  inaccurate  results 
in  such  cases  (cf.,  e.g.,  reference  3)*  The  skin  friction,  velocity 
profiles,  and  critical  Reynolds  number  are  calculated  for  several  Mach 
numbers  in  this  case  by  means  of  the  general  solution  derived  here. 

These  results  are  then  compared  with  known  series  solutions  of  the 
partial  differential  equations  (references  4  and  5)  for  this  case;  the 
agreement  is  shown  to  be  satisfactory  for  practical  purposes.  This 
agreement  includes  critical  Reynolds  numbers  for  laminar-flow  stability 
and  location  of  the  separation  point.  In  connection  with  the  latter  it 
is  significant  to  note  that  for  zero  Mach  number  the  Karman-Pohlhausen 
method  with  fourth-degree  velocity  profiles  leads  to  an  error  of  about 
30  percent  in  predicting  the  separation  point,  while  the  use  of  a  sixth- 
degree  profile  reduces  this  error  to  15  percent,  and  the  special  use  of 
a  seventh-degree  profile  for  this  (and  only  this)  purpose  practically 
eliminates  the  error  for  all  Mach  numbers  (O  to  3)  considered  numeri¬ 
cally  here. 

The  second  example  treated  is  the  flow  near  a  forward  stagnation 
point.  This  case  is  of  interest  not  only  because  it  may  represent  the 
subsonic  flow  over  a  blunt  nose  but  also  because  it  had  already  been 
found  by  Schlichting  and  Ulrich  (reference  6)  that  in  such  a  case  no 
physically  significant  values  can  be  found  for  the  boundary- layer  thick¬ 
ness  when  sixth-degree  profiles  are  used.  It  is  shown  in  the  present 
investigation  that  although  this  is,  strictly  speaking,  true,  an  approxi¬ 
mate  value  for  the  thickness  is  nevertheless  obtainable.  It  is  also 
shown  that  in  this  case  the  usual  type  of  application  of  the  Karman- 
Pohlhausen  method  with  fourth-degree  velocity  profiles  leads  not  only 


to  real  values  of  the  boundary-layer  thickness  but  also  to  results  of 
satisfactory  accuracy^  even  for  the  critical  Reynolds  number.  Conse¬ 
quently^  in  such  a  case  fourth-degree  profiles  may  be  used  with  confi¬ 
dence  in  the  accuracy  of  the  results  to  be  obtained. 

As  in  reference  1^  it  has  been  assumed  that  the  viscosity  coeffi¬ 
cient  is  directly  proportional  to  the  absolute  temperature^  and  that^ 
following  a  suggestion  of  Chapman  (references  7  and  8)^  the  constant  of 
proportionality  may  be  chosen  so  that  Sutherland’s  relation  is  exactly 
satisfied  at  the  wall.  This  introduces  considerable  mathematical  sim¬ 
plification  and  will  lead  in  practice  to  fairly  accurate  results  at 
least  for  lower  Mach  numbers  (below  5)-  Fon  the  same  reasons  it  has 
also  been  assumed  that  the  Prandtl  number  is  unity,  and  this  may  be 
viewed  as  an  approximation  for  air,  where  the  Prandtl  number  is  more 
nearly  0.72. 

This  work,  carried  out  at  the  Polytechnic  Institute  of  Brooklyn 
Aeronautical  Laboratories,  was  sponsored  by  and  conducted  with  the 
financial  assistance  of  the  National  Advisory  Committee  for  Aeronautics. 
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coefficient  of  t  in  velocity  profile;  also  given 
by  equation  (20) 

constant  "average"  value  of  a2 

constant  in  given  velocity  distribution  ^ 

equations  (27a)  and  (30a) 

factor  of  proportionality  in  equation  =  c(t/Tqq) 

skin-friction  coefficient  ^(p  Su/Sy)^  ^^p^u^^y2) j 

specific  heats  at  constant  pressure  and  at  constant 
volume ,  respectively 

integrals  defined  by  equations  (12) 

* 

values  of  F2,  and  F^  based  on  seventh-degree 

profiles;  given  by  equations  (B2) 
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constant  defining  strength  of  a  shock  wave  in  flow 
outside  of  boundary  layer;  given  by  equation  (l9e) 

integrating  factor  of  equation  (21) 

constant  defined  by  equations  (l8a)  and  (l8b) 

thermal  conductivity  of  fluid 

characteristic  length;  see  equations  (l3) 

characteristic  length  in  stability  calculations 
(appendix  C) 

Mach  number 

pressure 

Reynolds  number 

free-stream  Reynolds  number  ^p^u^L^p^) 

Sutherland’s  constant^  defined  by  equation  (8);  for 
air,  S  =  216°  R 

absolute  temperature 

transformation  variable  defined  by  equation  (lO) 

velocity  components  in  x-  and  y-directions ,  respectively 

coordinates  parallel  and  perpendicular  to  surface, 
respectively 

quantities  defined  by  equations  (C2) 
ratio  of  specific  heats  (^p/^v)^  alv,  7  =  1.4 
physical  boundary- layer  thickness  in  xy-plane 
boundary- layer  thickness  in  xt-plane 
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V 

I 

p 

T 

Subscripts 

o 

1 


value  of  X  based  on  seventh-degree  profiles  and  used 
only  for  determining  separation  point 

coefficient  of  viscosity 

kinematic  viscosity  (p/p) 

dimensionless  distance  along  wall  (x/^) 

mass  density 

dimensionless  variable  (t/S-^;;] 


values  at  wall,  for  example,  Tq 

local  values  at  outer  edge  of  boundary  layer;  for 
example,  M]_,  Tj,  and  u^ 


a 

b 

c 

s 


initial  value  in  region  of  adverse  pressure  gradient 
(appendix  B) 

values  immediately  behind  a  shock  wave  at  leading  edge 
values  at  critical  point  for  stability  (appendix  C) 
value  based  on  seventh-degree  profiles 


00 


values  in  undisturbed  free  stream;  for  example, 
T^,  and  u^ 


A  prime  ( ’ )  denotes  at  first  differentiation  with  respect  to  x. 
In  equation  (l4)  and  thereafter  in  the  main  text,  a  prime  denotes  dif¬ 
ferentiation  with  respect  to  In  appendix  C,  it  denotes  differen¬ 

tiation  with  respect  to  y/Z. 


BASIC  EQUATIONS 

Two-Dimensional  Compressible  Flow 


The  following  equations  characterize  the  steady  two-dimensional 
compressible  flow  of  a.  gas  in  the  laminar  boundary  layer  over  a  STirface 
of  large  radius  of  curvature  compared  with  the  boundary- layer  thickness: 

^  Su  .  ^  5u  ,  S  /  Su\ 


(1) 
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pucp 


Bx 


5t  • 


+  pvc 


P  ay 


-PlUiUi’u  + 


(2) 


pv) 

oy 


=  0 


(3) 


_P_ 

Pi 


Equations  (l)  and  (2)  are  the  momentui'i  and  energy  equations,  respectively. 
The  equation  of  continuity  is  given  hy  expression  (3)^  while  equation  (4) 
follows  from  the  ideal  gas  law.  In  equations  (l),  (2),  and  (4)  account 
has  been  taken  of  the  well-known  implication  of  Prandtl's  basic  boundary- 
layer  assumption  that  the  pressure  gradient  Bp/^y  across  the  boundary- 
layer  thickness  is  zero.  Thus  the  axial  pressure  gradient  can  be  con¬ 
sidered  as  given  by  the  Bernoulli  equation 


^  = 
Bx 


-PlUiUi 


(5) 


where  the  subscript  1  refers  to  the  potential  flow  over  the  surface 
streamline  outside  of  the  boundary  layer.  This  potential  flow  is 
assumed  to  be  known. 

For  a  Prandtl  number  pcp/k  of  unity  and  zero  heat  transfer  at  the 
wall  (Sf/Sy  =  0  at  y  =  0)  it  can  be  shown,  by  eliminating 
from  equations  (l)  and  (2),  that 


u2 


+ 


CpT 


=  Constant 


(6) 


Equation  (6)  can  be  expressed  in  the  following  nondimens ional  form: 


T 


Tl 


1  + 


7-1 


2 


(6a) 
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Thus  the  energy  equation  (2)  can  for  this  case  he  replaced  by  equa¬ 
tion  (6a),  which  gives  the  temperature  explicitly  as  a  simple  function 
of  the  velocity.  This  case  of  unit  Prandtl  number  and  zero  heat  trans¬ 
fer  is  the  one  which  will  be  treated  in  the  present  analysis. 

It  will  be  assumed  here,  as  in  references  1,  7,  and  8,  that  the 
viscosity  coefficient  can  be  expressed  as  a  linear  function  of  the 
temperature  in  the  form 


(7) 


where  C  is  a  constant  determined  so  that  relation  (7)  will  give  the 
same  value  of  |a  at  the  wall  as  the  Sutherland  viscosity-temperature 
relation.  The  latter  relation  is  generally  assumed  to  be  accurate  at 
least  below  the  hypersonic  range.  Thus 

^  ^  1  ^  (S;T„)  jg 

(Tq/Tco)  +  (s/tJ 

where^  S  is  a  constant  which  for  air  has  the  value  S  =  2l6°  R.  From 
equation  (6a)  it  follows  that  the  wall  temperature  Tq  as  a  function 
of  the  free-stream  Mach  nimiher  M^o  will  he: 


Hence  the  constant  C  will  he  a  function  of  the  Mach  number  M^  and 
the  free -stream  temperature 


By  integrating  equation  (l)  with  respect  to  y  across  the  boundary 
layer  thickness  y  =  0  to  y  =  6  and  introducing  the  Dorodnitzyn  vari¬ 
able  t  defined  by  the  transformation 


(10) 
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the  following  equation,  essentially  the  Karman  integral  moment lom  equa¬ 
tion  for  compressible  laminar  boundary  layers,  can  be  derived  (cf*  ref¬ 
erence  1,  equation  (24)): 


and  5-t^,  the  value  of  t  when  y  =  6,  is  the  boundary- layer  thickness 
in  the  xt -plane,  while  t  =  t/bf 

Introducing  the  dimensionless  variables  X  and  |  defined  by 


(13) 


where  L  is  a  characteristic  length  and  =  Pqo^ooVi^oo  1®  free- 
stream  Reynolds  number,  equation  (ll)  can  be  written  as 


X' 


+  2X 


Pi 


+ 


+ 


^1 


^  £2  h  ^ 

^1  Pi  ^00  ^1 


where  the  primes  henceforth  denote  differentiation  with  respect  to  | 
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Velocity  Profiles 

If  a  definite  form  of  profile  for  as  a  function  of  t  is 

assumed,  then  equation  (l4)  becomes  an  ordinary  differential  equation 
in  j.  In  the  present  analysis^  sixth-degree  profiles  will  he 
assumed  for  most  cases  and  will  he  chosen  to  satisfy  appropriate  bound¬ 
ary  conditions  at  the  wall  and  at  the  outer  edge  of  the  boundary  layer. 
This  is  in  accordance  with  the  results  of  reference  1,  where  it  was 
concluded  that^  with  the  possible  exceptions  of  determination  of  sep¬ 
aration  point  in  an  adverse  pressure  gradient  and  of  calculation  of  the 
boundary  layer  in  the  vicinity  of  a  stagnation  pointy  the  use  of  sixth- 
degree  profiles  should  for  practical  purposes  lead  to  sufficiently  accu¬ 
rate  deteimiination  of  laminar-boundary-layer  characteristics,  including 
critical  Reynolds  number  for  stability.  Sixth-degree  profiles  have  also 
been  applied  recently  by  Weil  (reference  9)- 

The  following  boundary  conditions,  expressed  in  the  xt-plane,  will 
be  satisfied: 


At  T  =  0: 


At  T  =  1: 


=  0 

^1 

It-i  /  p  \u  ]  T 

^  Q 

=  1 

^1 

St  ^^3 


(15) 


It  should  he  noted  that  the  second  condition  at  the  wall  (t  =  0)  in  equa¬ 
tions  (15)  follows  in  general  from  the  original  partial  differential  equa¬ 
tion  (1)  in  conjunction  with  equation  (7).  The  third  condition  at  the 
wall  follows  from  differentiation  of  equation  (l)  with  respect  to  t,  in 

conjunction  with  the  condition  of  zero  heat  transfer  ((St/St)q  =  o)  at  the 
wall.  The  second  and  third  conditions  of  equations  (15)^  moreover,  are 
valid  only  for  zero  normal  velocity  at  the  wall. 
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The  sixth-degree  polynomial  satisfying  conditions  (15)  is  found 

to  be 


—  =  2(1 

ui  \ 


^  2 
— jr  +  agT 


-  ^5  +  2a2)T^  +  2^3  +  82) 


where: 


a2 


2  ^co  P„\Ti/ 


(17) 


Substitution  of  equation  (I6)  into  equations  ( 12)  leads  to  the 
following  explicit  expressions  for  F2^  and  F^  in  terms  of  82: 


Fj^  =  0.1093  +  0.0021132  -  0. 00062282^ 

F2  =  -0.2857  -  0.0190532  -  (0.3950  + 

0. 0211632  -  0. 00062282^) 


(18) 


F^  =  2  -  0. 400082 


By  substituting  expressions  (18)  and  (l7)  into  equation  (l4);  the  latter 
becomes  an  ordinary  differential  equation  of  the  first  order  in  X(|). 

The  quantities  with  subscript  1  may  be  regarded  as  given  functions 
of  I  from  the  potential  flow  outside  of  the  boundary  layer.  From  the 

relation  (u2/2)  +  CpT  =  Constant,  which  is  valid  outside  of  the  boundary 
layer  even  across  a  shock  wave,  it  follows  that 


(19a) 
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Moreover,  the  local  Mach  number  M3_  can  he  expressed  by; 


Ml 


2 


(I9h) 


If  there  are  no  shock  waves ^  the  flow  outside  of  the  boundary  layer  may 
he  considered  as  isentropic^  so  that  the  relation 


Pi  _  /TiU/Cr-i) 

P  \T  / 

\  00  / 


(19c) 


Is  valid.  In  the  presence  of  a  shock  wave,  the  flow  may  be  considered 
as  isentropic  before  and  behind,  though  not  across,  the  wave.  In  that 
case,  equation  (l9c)  may  be  modified  to 


Pi  /Ti\i/(7-1) 

Poo  “  vJ 


(I9<i) 


where  G  is  a  constant  given  by 


G  = 


yPeoJ  Y^oo  / 


(I9e) 


while  pi3  and  T-j^  are  the  values  of  p  and  T  immediately  behind  the 
wave^  and  *  *^1  "to  the  region  behind  the  wave  at  the  outer 

edge  of  the  boundary  layer.  Relation  (l9c)  can  be  considered  as  a  special 
case  of  relation  (l9d);  namely  the  case  of  G  =  !•  The  ratios  P^/Poo 

T-j^^Too  may  be  considered  as  known  from  the  given  flow  outside  of  the  bound 
apy  layer. 

By  the  use  of  equations  (9),  (l9a)^  and  (l9d.)  the  expression  (l7)  for 
32  can  be  written^  more  explicitly^  as 


3-27 
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GENERAL  SOLUTIONS  OF  EQUATIONS 


For  a  given  potential-flow  velocity  distribution  u^/uj  0  and 
Mach  number  M^>  the  differential  equation  (l4)  can  be  solved  for  X( I) 

without  any  basic  difficulties  by  well-known  numerical  methods,  such 
as  Adam's  method  or  the  Runge-Kutta  method  (cf.,  e.g.,  reference  10). 
Such  a  straightforward  procedure  may  in  practice  nevertheless  be  tedious. 
A  relatively  simple  general  approximate  solution  of  equation  (14)  will 
therefore  be  derived  in  this  section.  This  solution  will  be  found  to 
be  sufficiently  accurate  for  most  practical  purposes. 


General  Approximate  Solution 

The  solution  to  be  derived  here  will  be  based  on  the  approximating 
ass\imption  that  in  the  expressions  for  F]_  and  F2  the  quantities  a2 
and  a2^  as  given  by  equation  (20)  may  be  replaced  by  constant  average 

values  ig  and  ag^  for  the  flow.  This  assumption  is  justified  by  the 
fact  that  in  expressions  (18)  for  Fj_  and  Fg  the  ag  terms  are  rela¬ 
tively  small;  hence  even  large  errors  in  the  evaluation  of  ag  will  lead 
to  only  small  errors  in  the  evaluation  of  F^  and  Fg. 

By  replacing  ag  and  ag^  by  constant  average  values  ag  and 
in  expressions  (I8)  for  Fi  and  Fg,  the  quantity  F^  becomes  a 
constant,  while  Fg  can  be  written  as 


Fg  =  -0.2857  -  0.01905ag 


(l8a) 


where  K  is  a  constant  given  by 


K  =  0.3950  +  0. 02116a g 


0. 000622a 


(l8b) 
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By  use  of  the  relation 


which  follows  from  equation  (6a),  and  hy  insertion  of  expression  (l8a) 
for  F2^  expression  (l8)  for  and  expression  (l?)  for  82  in  the 

right  side  of  equation  (l4),  the  latter  equation  can  he  written  in  the 
form: 


+  -^(0.08571  +  0.0190582)  + 

Ff  \  * 


Jl  ^  ^ 

Fi  ui  Pi  T„ 


(21) 


Equation  (21)  is  now  a  linear  differential  equation  of  the  first 
order  in  A,.  By  use  of  relations  (l9a),  (l9'b),  and  (l9d)  it  will  be 
found  that  an  integrating  factor  I  of  equation  (2l)  is: 


I 


^  (0.08571+0 . 0190552)/  T 


The  solution  of  equation  (21)  satisfying  the  boundary  condition  X  =  0 
or  some  finite  value  (in  the  case  of  a  stagnation  point  at  |  =  O)  at 
1=0  is  then: 


^  Jo 


^  ^ (0*0^571+0.0190532)^ i(o. 1950+0. 0211 6a2-0,000622l^  ) 


Vic 


dl 


FiG 


,  +  ^(o* 08571+0. 01905a2L\-^-  ^(0.1950+0.0211632-0.000622^) 

fey 


(22) 


where  T]_^T^  is  obtained  from  equation  (l9a). 
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For  a  given  velocity  distribution  effect  of  Mach 

number  is  given  by  the  terms  in  equation  (22)  and^  in  the 

case  of  a  shock  wave  at  the  leading  edge,  also  by  the  constant  G 
(cf.  equation  (l9e)).  If  the  flow  over  a  given  object  is  considered, 
however,  then  the  Mach  number  effect  will  be  contained  also  in  the 
Ui|Uoo  terms,  since  the  velocity  distribution  will  then  be  a 

function  of  Mach  nxxmber. 

The  integral  in  equation  (22)  can  in  practice  be  evaluated  without 
difficulty  by  niomerical  means  such  as  Simpson’s  rule  (cf.,  for  example, 
reference  lO).  An  average  value  "32  for  B2  can,  in  any  given  case, 

be  chosen  by  considering  equation  (20). 


Boundary-Layer  Characteristics 

Once  \( I)  has  been  determined,  the  boundary -layer  characteristics 
follow  from  the  equations  developed  here.  Thus  the  local  skin-friction 
coefficient  Cf  will  be;! 


Cf 


(>.au/ay)o  Ti  (^  -  t)  ^  ,i/3 

p  u  ^/2 

^00  00  / 


(23) 


The  physical  "boundary-layer  thickness  in  the  xy-plane  can  he  determined 
from  transformation  (lO),  and  an  explicit  expression  for  it  in  terms  of 
X  is  given  in  appendix  A. 

The  velocity  profiles  in  the  xt-plane  are  given  hy  equation  (l6). 
These  can  he  converted  to  the  physical  xy-plane  hy  means  of  transforma¬ 
tion  (10).  An  explicit  expression  for  y  in  terms  of  t  is  given  in 
appendix  A.  The  temperature  profiles  follow  from  the  velocity  profiles 
hy  means  of  equation  (6a).  It  should  he  noted,  however,  that  for  most 
practical  purposes  it  is  not  actually  necessary  to  determine  the  profiles 
in  the  xy-plane.  Skin  friction  (cf.  equation  (23))  and  separation  point, 


^It  is  demonstrated  subsequently  that  the  flow  near  a  stagnation 
point  is  better  represented  by  a  fourth-degree  velocity  profile  than  hy 
a  sixth-degree  profile.  In  such  a  region  the  expression  for  the  skin- 

I  ®2 

friction  coefficient  is  altered  in  that  the  term  11  - 

f.  ®2\ 


L  IC 

1 


is  replaced 
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for  example^  can  “be  determined  directly  from  results  in  the  xt-plane^ 
while  stability  calculations  (cf.  appendix  C)  can  be  performed  without 
use  of  the  profiles  in  the  xy-plane. 

The  equations  developed  thus  far  should  be  adequate  for  the  prac¬ 
tical  determination  of  boundary- layer  characteristics  in  most  cases* 
There  are,  however,  two  important  exceptions.  These  are:  (a)  The  loca¬ 
tion  of  the  separation  point  in  an  adverse  pressure  gradient  and  (b)  the 
flow  in  the  vicinity  of  a  forward  stagnation  point.  Item  (a)  is  treated 
in  the  immediately  succeeding  paragraph,  while  item  (b)  is  investigated 
in  a  later  separate  section. 


Determination  of  Separation  Point 

Although  the  equations  thus  far  developed  will  usually  be  found  to 
lead  to  sufficiently  accurate  results  for  most  boundary- layer  character¬ 
istics  even  in  flows  with  adverse  pressure  gradients,  it  will  be  found 
that  the  location  of  separation  points  will  still  be  predicted  with 
appreciable  error.  This  error  can  be  greatly  diminished  by  a  simple 
modification  of  the  preceding  equations. 

For  the  purpose  of  determining  the  separation  point,  velocity  pro¬ 
files  can  be  chosen  which,  in  addition  to  satisfying  boundary  condi¬ 
tions,  (  15)^  satisfy  a  further  boundary  condition  at  the  separation  point. 
The  additional  boundary  condition  is  obtained  by  differentiating  equa¬ 
tion  (l)  twice  with  respect  to  t,  and  by  then  taking  values  at  the  wall 
at  the  point  where  Su/St  =  0.  In  this  manner,  with  the  use  of  equa¬ 
tion  (6a),  the  condition 

S^(u/uij 


is  obtained. 

Condition  (24)  is  strictly  valid  only  at  the  separation  point  (where 
Su/5y  =  O).  Nevertheless,  by  using  this  condition  for  the  entire  flow, 
it  seems  plausible  that  the  location  of  the  separation  point  in  any  given 
case  could  be  predicted  with  improved  accuracy,  since  the  value  of  the 
dimensionless  boundary- layer  thickness  X  may  be  expected  to  be  thereby 
more  accurately  determined  at  the  separation  point.  This  has  already 
been  shown  by  an  example  for  incompressible  flow  in  reference  2  and  will 
be  further  shown  by  examples  (table  I)  for  compressible  flow  in  the  present 
analysis. 


=  0 


T=0 


(24) 
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A  seventh-degree  polynomial  in  t  satisfying  conditions  (15)  and 
(214-)  can  he  chosen.  By  proceeding  in  the  same  manner  as  (previously) 
with  the  sixth-degree  profiles,  the  following  solution  X,g  for  X(5), 

analogous  to  equation  (22),  is  obtained  (a  detailed  derivation  is  given 
in  appendix  B) ; 


= 


32.83  Jo 


/^\6. 128/^X0.9366 

\  \  “  / 


/u^\7.123/^\2.437 

j  yToo  j 


(25) 


The  value  of  y  =  1.4  has  been  used  in  equation  (25).  In  case  the 
region  of  adverse  pressure  gradient  follows  a  region  of  favorable  pres¬ 
sure  gradient,  and  hence  starts  at  some  point  behind  the  leading  eige, 
then  equation  (25)  should  be  modified  to  equation  (b4)  in  appendix  B. 
Separation  will  occur  at  the  point  where  (see  appendix  B): 


The  location  of  the  separation  point  is  thus  determined  at  the 
value  of  I  for  which  equations  (25)  (or  (b4))  and  (26)  give  the 
same  value  of  Jig.  For  any  given  adverse  pressure  gradient  (negative 

Ui'({))  and  free-stream  Mach  number  the  function 

plotted  against  |  in  accordance  with  both  equation  (25)  (or  equa¬ 
tion  (b4))  and  equation  (26)  in  the  (anticipated)  vicinity  of  the 
separation  point.  The  point  of  intersection  of  these  two  curves  is 
then  the  separation  point It  should  he  emphasized  that  equation  (25) 


%his  point  will,  according  to  equations  (25)  and  (26),  be  inde¬ 
pendent  of  G;  hence  one  may  put  G  -  1  for  this  purpose.  This  does 
not  imply  that  the  separation  point  will  actually  be  unaffected  by  a 
shock  wave  at  the  leading  edge,  since  the  values  of 

U2_’/u^(|)  may  still  contain  the  effect  of  such  a  wave. 
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is  to  be  used  only  for  determination  of  separation  point.  For  other 
purposes,  such  as  determination  of  .skin-friction  distribution  or  of 
critical  Reynolds  number,  equation  (22)  and  its  related  equations  are 
to  be  used,  even  for  a  region  in  an  adverse  pressure  gradient,  unless 
properties  close  to  or  at  the  separation  point  are  sought. 

An  application  of  the  general  solutions  developed  here,  as  well  as 
a  check  on  the  accuracy  of  results  obtained  from  them,  will  now  be  given 
by  means  of  examples  based  on  two  different  types  of  flow  outside  of  the 
boundary  layer. 


FLOW  WITH  LINEARLY  DIMINISHING  VELOCITY 


The  case 


(27a) 


where  b  is  a  positive  constant,  is  now  considered.  This  represents 
the  simplest  type  of  adverse  pressure  gradient.  By  introducing  the 
linear  change  of  variables 


J 


(28) 


equation  (l4)  remains  unchanged,  except  that  i  is  replaced  by 
and  X,  by  Moreover,  ■a]_/u^  now  becomes 


^1 


u. 


=  1-1- 


(2Tb) 


Thus  the  constant  b  no  longer  appears  in  the  equations.  Consequently, 
in  this  case  it  is  actually  necessary  to  solve  equation  (l4)  only  for 
b  =  1,  since  for  any  other  value  of  b  it  is  necessary  only  to  replace  | 
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ty  li  and  \  by  in  the  solution  obtained  for  b  =  1.  The 

present  example^  therefore,  will  treat  the  case 

^  =  1  -  I  (2Tc) 

^co 

A  constant  average  value  32  '^or  82  can  be  chosen  by  noting  that 
in  this  case  82  will  vary  from  0  (at  |  =  O)  to  5  (at  the  separation 
pointy  where,  according  to  equation  (16),  1  -  (^2/^)  ”  Consequently 

the  value  ^2  =  ^*5  niay  be  chosen.  Moreover,  the  value  a^  =  12.5  will 

be  used  here.  With  these  values,  the  approximate  solution  (22)  for  the 
case  represented  by  equation  (27c)  without  a  shock  wave  (G  =  l)  at  the 
leading  edge  becomes: 


(29a) 


where  t]  =  1  -  For  zero  Mach  number,  that  is,  =  0,  equation  (29a) 
reduces  to: 


X  =  5.76  (1  -  3 


(29b) 


The  skin  friction  follows  from  equations  (23)  and  (20),  while  the 
velocity  profiles  can  be  obtained,  as  previously  explained,  by  means 
of  equations  (16)  and  (A2)  (appendix  A).  For  M^  =  1  and  3^  the  inte¬ 
gration  in  equation  (29a)  was  performed  by  expanding  the  integrand  in 
a  series  according  to  the  binomial  theorem.  The  results  for  X.(^), 
skin  friction  at  the  wall,  and  velocity  profiles  are  shown  in  figures  1, 
2,  and  3^  respectively,  with  M^^  =  0,  1,  and  3. 
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The  critical  Reynolds  number  for  the  stability  of  the  laminar 
boundary  layer  in  the  presence  of  small  disturbances  was  calculated 
here  on  the  basis  of  criteria  developed  by  Lin  and  lees  (references  11 
and  12).  The  method  of  calculation  is  outlined  in  appendix  C,  and  the 
results  for  the  critical  Reynolds  number  PoqU^l/p^  at  the  station 

I  =  0.o496  are  given  in  table  II  for  =  0,  1,  and  3.  By  comparing 
the  results  for  zero  Mach  number  with  those  for  the  critical  Reynolds 
number  (1,500,000)  at  the  same  station  for  the  flow  over  a  flat  plate 
based  on  the  solutions  of  reference  8  (see  also  reference  l),  it  is 
seen  that  the  adverse  pressure  gradient  is  here  highly  destabilizing. 
Moreover,  it  can  also  be  seen  from  table  II  that  higher  Mach  numbers 
greatly  destabilize  the  flow  over  an  adiabatic  wall. 

As  a  check  on  the  accuracy  of  approximate  solution  (22)  to  differ¬ 
ential  equation  (l4),  equation  (l4)  was  solved  exactly  for  this  case 
for  M^  =  0  and  numerically  by  Adam's  parabolic  method  (reference  13) 
with  increments  of  A|  =  0.01  for  /  0.  The  solutions  thus  obtained 

for  X(5)  and  for  the  profiles  are  shown  in  figures  1  and  3,  respec¬ 
tively.  The  agreement  between  these  solutions  and  those  based  on  equa¬ 
tion  (22)  is  seen  from  these  figures  to  be  satisfactory  for  practical 
purposes.  As  a  further  check,  the  critical  Reynolds  number  at  |  =  0.0496 
was  calculated  for  Moo  =  0  by  using  the  exact  solution  of  equation  (l4). 
Table  II  indicates  satisfactory  agreement  between  this  result 
((®oo)c  380,000)  and  that  ((Roo)c  "  330, OOo)  based  on  equation  (22). 

It  is  thus  indicated  that  equation  (22)  is  a  satisfactory  approxi¬ 
mation  to  an  exact  solution  of  differential  equation  (l4).  It  is  also 
desirable  to  check  whether  an  exact  solution  of  equation  (l4),  or 
approximation  (22),  is  a  satisfactory  approximation  to  an  exact  solu¬ 
tion  of  the  original  partial  differential  equation  (l).  This  would  con¬ 
stitute  a  check  on  the  practical  reliability  of  the  basic  method  used 
here,  namely  the  extension  of  the  Karmn-Pohlhausen  method  to  sixth- 
degree  profiles.  Such  a  check  can  be 'made  by  comparing  the  results 
obtained  by  means  of  equations  (l4)  and  (22)  with  the  results  of  Howarth 
(reference  4)  for  zero  Mach  number,  based  on  an  accurate  solution  of  the 
original  partial  differential  equations  (l)  and  (3)»  Such  a  comparison 
is  shown  in  figures  2  and  3  for  skin  friction  and  velocity  profiles, 
respectively,  and  in  table  II  for  critical  Reynolds  niimber.  In  all 
cases  the  comparison  indicates  that  the  solutions  obtained  by  means  of 
the  equations  developed  here  give  sufficiently  reliable  results  for 
practical  p\rrposes. 

The  separation  point  is  determined  by  equating  the  right  sides  of 
equations  (25)  and  (26).  In  the  present  case  the  point  at  which  these 
right  sides  are  equal  was  determined  by  plotting  Xg(5)  according  to 
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both  equations  (25)  and  (26)  on  the  same  sheet  in  the  vicinity  of 
^  =  0.l4  for  Hx)  =  0  and  1  and  |  =0.10  for  =  3.  These  are  the 

vicinities  in  which,  according  to  figure  1  (based  on  sixth-degree,  not 
seventh-degree,  profiles),  separation  might  be  expected  to  occur.  The 
separation  points  determined  by  this  means  for  =  0,  1,  and  3  are 
shown  in  table  I,  together  with  the  results  of  Stewartson  (reference  5). 
The  agreement  between  these  two  sets  of  results  is  seen  to  be  excep¬ 
tionally  good. 

For  purposes  of  comparison,  the  separation  point  was  also  deter¬ 
mined  by  the  use  of  equation  (22)  and  by  means  of  the  classical  Kanrian- 
Pohlhausen  method  based  on  fourth-degree  profiles,  the  latter  for  zero 
Mach  number  only.  The  results  are  shown  in  table  I.  Although  the  use 
of  the  sixth-degree  profiles  (equation  (22))  gives  a  more  accurate 
location  of  the  separation  point  than  the  fourth-degree  profiles,  the 
special  use  of  the  seventh  degree  (equations  (25)  and  (26))  for  this 
purpose  gives  here  virtually  the  exact  values. 

It  may  be  asked  whether  the  seventh-,  instead  of  the  sixth-degree 
velocity  profile  could  be  used  for  determining  the  other  boundary -layer 
characteristics,  in  addition  to  the  separation  point.  Consequently,  the 
skin  friction  and  velocity  profiles  for  the  present  case  were  determined 
on  the  basis  of  the  seventh- degree  velocity  profile,  in  conjunction  with 
equation  (l4)  (cf.  appendix  B).  It  was  found  that  the  results,  even  at 
a  station  fairly  close  to  the  separation  point,  did  not  agree  so  closely 
with  Howarth*s  solution  (reference  4)  as  those  based  on  the  sixth-degree 
profile.  It  is  therefore  concluded  that  the  seventh- degree  profile 
should  be  used  only  as  a  means  of  determining  the  separation  point  in 
regions  of  adverse  pressure  gradients  and  that  otherwise  (except  for 
flow  near  a  stagnation  point)  the  equations  developed  on  the  basis  of 
the  sixth-degree  velocity  profile  should  be  used. 


FLOW  MIAR  A  FORWARD  STAGNATION  POINT 


The  case 


u. 


(30a) 


represents  physically  the  flow  in  the  vicinity  of  a  forward  stagnation 
point,  such  as  over  the  leading  edge  of  a  hlunt  object  in  subsonic  flow. 
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As  in  the  preceding  case  (equation  (27a))  the  constant  h  here  can  he 
eliminated  hy  means  of  transformations  (28).  Hence  only  the  case 


u 


00 


(30b) 


need  be  treated. 

For  zero  Mach  number  =  o)  an  exact  solution  of  ordinary  differ¬ 
ential  equation  (l4)  in  the  case  represented  by  equation  (30b)  is  of  the 
form  X  =  Constant.  This  can  also  be  seen  from  solution  (22)  for  this 

case,  with  ^2^  ~  where  according  to  equation  (20) 

(with  Moo  =  0,  G  =  1,  and  u^'/u^  ~  ^2  “  From  equation  (22), ^ 

in  fact,  the  following  cubic  in  X  is  obtained; 


f(x)  =  X^  +  37.434^  -  979. 4x.  +  6435  =  0 


(31) 


By  plotting  f(X.)  against  X,  it  is  fotmd  that  there  are  no  physi¬ 
cally  significant  roots  of  equation  (31).  Such  a  result  was  already 
obtained  by  Schlichting  and  Ulrich  (reference  6).  It  will  be  found, 
however,  that  the  ciirve  f(X)  against  X  in  the  physically  significant 
region  comes  relatively  close  to  the  \-axis  at  a  point  when  the  curve 
has  a  local  maximum  value.  If  the  value  of  A,  at  this  point  is  taken 
for  purposes  of  an  approximation,  then 


X  =  9.481 


(32) 


Since  the  approximation  (32)  does  not  appear  to  be  satisfying  in 
principle,  the  use  of  fourth-degree  velocity  profiles  will  be  inves¬ 
tigated  for  this  case.  The  fourth-degree  polynomial  satisfying  all  the 
boundary  conditions  in  equations  (15)  with  the  exception  of  the  two  con¬ 
ditions  involving  the  third  derivative  of  velocity  j =  0 

at  T  =  0  and  at  t  =  1^  is 
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With  the  profile  (33) ^  equation  (l4)  remains  valid^  but  the  F's  as 
defined  by  equations  (12)  are  now  given  explicitly  by  the  following 
expressions: 


=  0.1175  +  0.00212a2  -  O.OOOl+^lag^ 

Fg  =  -0.3000  -  0.0166732  -  ^  ~  -  M3^2^o.i*-175  + 
O.Ol878a2  -  0.00041Ha2^)  [ 


(3^) 


F3  =  2  -  0.333382 


Expressions  (34)  replace,  in  this  case,  expressions  (18)  based  on  sixth- 
degree  profiles.  For  the  case  (30b)  with  =  0,  equation  (l4)  leads 

to  a  cubic  in  X  whose  physically  significant  root  is  now  found  to  be 
(cf.  also  reference  3) 


X  =  7.052 


(35) 


To  compare  the  accuracy  of  solutidn  (32),  based  on  sixth-degree 
velocity  profiles,  with  that  of  solution  (35) ^  the  skin  friction,  veloc¬ 
ity  profiles,  and  critical  Reynolds  number  for  laminar-flow  stability 
have  been  calculated  on  the  basis  of  both  of  the  solutions.  The  results 
thus  obtained,  together  with  a  comparison  with  results  of  an  exact  solu¬ 
tion  (reference  l4)  for  this  case,  are  shown  in  tables  III  and  IV  and  in 
figure  4.  From  these  comparisons  it  is  seen  that  the  solution  based  on 
fourth-degree  profiles  leads  in  this  case  to  more  accurate  results  than 
that  based  on  sixth-degree  profiles  for  the  boundary- layer  character¬ 
istics  considered,  including  the  critical  Reynolds  number.  It  may  there¬ 
fore  be  concluded  that  for  a  stagnation-flow  case  represented  by  equa¬ 
tions  (30a)  and  (30b),  a  satisfactory  solution  for  practical  purposes 
can  be  obtained  by  means  of  equations  (33) >  (34),  and  (35)  based  on 
fourth-degree  velocity  profiles. 
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For  ^  0,  an  approximate  solution  of  equation  (l4)  quite  similar 

to  equation  (22)  can  he  obtained  for  the  present  case.  Using  expres¬ 
sions  (3^)  and  putting  =  I  in  accordance  with  equation  (3013)^ 

the  following  approximate  solution  is  thus  obtained:  - 


F^G 


\  = 


Uo 


3+  f^(o*1333+O.Ol667a2)^rj,^\-^-  ^(0.2508+0. 0187852-0. 00044la2^) 


4+  ^(0. 1333+0. 0l667a2).^^.-^  - ^(0.2508+0.0187852-0. 00044182^) 


(36) 


where  is  now  given  by  equation  (3^).  At  the  leading  edge  (|  ==  O) 

the  value  Xq  of  X  will  be,  according  to  equation  (36), 


7-2 

^7-1 


(0.1333  +  0.0l667a2  + 


(36a) 


As  a  numerical  example,  the  case  =  1  was  calculated  by  means 

of  equations  (36)  and  (36a),  assuming  a’2  =  -4.48  and  3,2^  =  20.1.  The 

results  for  X(|),  skin  friction,  and  velocity  profiles  are  shown  in  fig¬ 
ures  5^  6,  and  7>  together  with  the  results  for  =  0. 


CONCLUSIONS 


From  the  foregoing  analysis,  based  on  an  extension  of  the  Karm^n- 
Pohlhausen  method  primarily  to  sixth-degree  velocity  profiles  for  deter¬ 
mining  laminar-boundary-layer  characteristics  in  the  compressible  flow 
over  an  adiabatic  wall  in  the  presence  of  an  axial  pressure  gradient, 
the  following  main  conclusions  can  be  drawn. 

1.  For  a  given  free-stream  Mach  number  and  a  given  velocity  distri¬ 
bution  outside  of  the  boundary  layer,  a  relatively  simple  general  approxi¬ 
mate  solution  to  the  boundary- layer  equations  has  been  derived  for  a 
Prandtl  number  of  unity  and  for  a  linear  relation  between  the  viscosity 
coefficient  and  the  temperature.  This  solution  not  only  is  convenient 
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to  apply  but  also  leads  to  results  which  will  be  sufficiently  accurate 
for  most  practical  purposes,  including  determination  of  critical  Reynolds 
number. 

2.  For  the  particular  purpose  of  determining  the  location  of  the 
separation  point  in  a  region  of  adverse  pressure  gradient,  a  relatively 
quick  and  simple  method  has  been  given  here,  based  on  the  use  of  seventh- 
degree  velocity  profiles.  The  numerical  example  indicated  that  the 
results  thus  obtained  should  in  general  be  quite  accurate. 

3.  In  the  special  case  of  flow  near  a  forward  stagnation  point,  the 
Karman-Pohlhausen  method  with  fourth-degree  velocity  profiles  gives 
results  of  very  satisfactory  accuracy,  even  for  the  critical  Reynolds 
number. 

4.  The  numerical  examples  given  here  illustrate  the  destabilizing 
influences  of  an  adverse  pressure  gradient  and  of  a  Mach  number  greater 
than  zero. 


Polytechnic  Institute  of  Brooklyn 
Brooklyn,  N.  Y.,  May  28,  1951 
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APPENDIX  A 

PHYSICAL  BOUNDARY-LAYER  THICKNESS 


Results  in  the  xt-plane  can  be  transformed  into  the  physical 
xy-plane  by  means  of  relation  (10),  Thus 


(Al) 


By  inserting  velocity  profile  (16)  and  equations  (l9a)  and  (l9b)  into 
expression  (6a)  for  equation  (Al)  yields: 


82(1  - 

t5  . 

-1(5. 

2a2)(l  -  ^)t6  . 

®2)(l 

a2\ 
5  / 

2 

“  7 
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+ 

?(3 

+  82)  - 

K- 

■  Tj 

82)]  P  + 

lE’" 

232)=^ 

-  282  ( 

f 

2  + 

< 

.9  _ 

82)(5 
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10 

P 

+  232) ( 

2.| 
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11(3 

^  11 
+  82]^  T  - 

i(3  .  a,)(.  .  I  .  I  3,)^«  - 


(A2) 
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After  has  been  obtained,  the  point  5,y  in  the  physical  plane 

corresponding  to  any  given  point  in  the  mathematical  plane  can  be 

directly  determined  by  means  of  equation  (A2) . 

The  boundary-layer  thickness  5  in  the  physical  plane  is  the  value 
of  y  at  which  t  =  1.  Hence,  by  putting  t  =  1  in  equation  (A2),  the 
following  explicit  expression  is  obtained  for  the  boundary-layer 
thickness:  . 


(7  -  1)  ^  2/Ui\2 


^  =  1  + 


w 


1  + 


1  - 


^1' 


u 


(0.3950  +  0. 0211682  -  0.00062282^)  (A3) 
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APPENDIX  B 

DETERMINATION  OF  SEPARATION  POINT 


A  seventh-degree -polynomial  velocity  profile  satisfying  condi¬ 
tion  (24)  in  addition  to  conditions  (15)  is 


“  (l  ■  ^  *  3=2)^^  -  (I  +  *2)^''  (Bi) 


where  is  given  hy  equa.tion  (20).  With  the  profile  (Bl),  the  F’s 

as  defined  hy  equations  (12)  hecome  (a  subscript  s  is  used  for  the 
present  results  to  distinguish  them  from  the  corresponding  results  based 
on  sixth-degree  profiles); 


F^g  =  0.1156  +  0.0025392  -  O.OOl454a2^ 

Fgg  =  -0.0297682  -  0.3125  -  ^^^2:  Mi2^0.4281  + 

0.0322992  -  0.0014532^) 


(B2) 


F^g  =  1.7500  -  0.5000a2 


The  separation  point  will  be  located  where 


it 


therefore  follows  from  equation  (Bl)  that  laminar-flow  separation  will 
occur  at  the  point  where  =  3«5*  By  using  expressions  (B2)  and 


assxaming^  as  in  the  section  entitled  ^'General  Approximate  Solutions^ 


V 
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that  in  the  expressions  for  F]_s  and  F2s  quantities  82  and 

may,  as  an  approximation,  be  replaced  by  constant  values  a'2 
the  following  approximate  solution  for  A.(0  of  differential  equa¬ 

tion  (l4)  (with  F-^^  F2^  and  F^  replaced  by  F^^g,  F2s^  and  F^g, 
respectively)  is  obtained: 


K  =  - 

s  2  F 


7  1  Jo  v^) 


Is^  +^{0. 0625+0 . 029T6a2)^Ti'^^  *“  (^ •  •  03229a2-0 .  OOl45a2^  ) 


(B3) 


Since  the  main  oojecu  01  equation  (B3)  is  to  determine  the  separa¬ 
tion  point,  the  values  of  82  and  may  for  this  purpose  be  chosen 

as  those  which  82  and  82^  would  actually  have  at  that  point.  Hence, 
for  determining  the  separation  point,  the  values  82  =  3.5  and 
82^  =  12.25  are  inserted  into  equation  (B3).  With  7  =  1.4,  expres¬ 
sion  -(25)  for  Xs(^)  then  results. 

Since  separation  occurs. where  82  =  3-5  and  since  82^  which  is 

proportional  to  X,  is  given  by  equation  (20),  the  value  of  X  at  which 
separation  occurs  will  be  that  given  by  equation  (26)  in  the  main  text. 

Equations  (B3)  and  (25)  for  Xg(^)  are  valid  only  under  the  bound¬ 
ary  condition  Xg  =  0  at  |  =  0,  In  case  the  region  of  adverse  pressure 
gradient  starts  at  some  point  ^  other  than  the  leading  edge,  this 

boundary  condition  must  be  replaced  by  the  condition  Xs  =  Xg  at  I  =  ^g. 

Solution  (25)  for  Xg(l)  must  then  be  replaced  by  the  following  solution 
fwith  a2  =  3.5  and  82^  =  12. 25V 
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The  value  of  is  obtained  directly  from  equation  (22)  of  the  main 

text •  Thus 


= 


>a  /  \3  + 


4  Jo  VH 


^(0.08571+0.0190532)/  -i(o.  1950+0. 0211632-0.0006223.2^) 


a  FnG 


^(0.085Tl+0.01905a2)^TiN^y§3--i(0.1950+0. 0211652-0. 


000622a. 


) 


(B5) 
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APPENDIX  C 


DETERMINATION  OF  CRITICAL  REYNOLDS  NUMBER 


In  accordance  with  the  procedure  outlined  hy  Lees  in  reference  11^ 
the  critical  Reynolds  number  at  a  given  station  x  for  the  stability 
of  the  laminar  boundary  layer  can  be  determined  in  the  following  manner. 
It  is  required  first  to  determine  the  value  of  y/Z  such  that 


where 


the  prime  here  denotes  differentiation  with  respect  to  y/Z  and  Z  is  a 
convenient  length.  In  the  present  analysis  the  quantity  6^  was  chosen 

for  Z^  except  in  the  case  of  Howarth's  series  solution  (reference  k)  of 

the  partial  differential  equations,  where  the  quantity  was 

chosen  for  Z. 

The  appropriate  value  of  y/ Z  satisfying  equation  (Cl)  must  usually 
be  found  graphically  by  trial.  After  this  value  has  been  determined,  the 
critical  Reynolds  number  based  on  Z  and  on  local  values  can  be  obtained 
from  the  expression: 
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where  the  subscript  c  denotes  values  at  the  point  where  equation  (Cl) 
is  satisfied- 

For  compressible  flow  it  is  more  convenient  to  perform  calculations 
in  the  xt-plane  than  in  the  xy-plane.  With  Z  =  5^^  the  foregoing  rela¬ 
tions  can  be  expressed  directly  in  the  xt-plane  by  use  of  the  relations: 

5  _  hr  /  S  \  _  ^1  h 

S(y/0\ST/  T  St 


(’hY  ii  ^ 

where  t  =  t/6^.  Using  a  dot  (*)  to  denote  differentiation  with  respect 
to  Ty  equations  (C2)  and  (C3)  become 
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The  quantity  y/6^  in  terms  of  t  is  given  hy  equation  (A2)  in 

appendix  A>  the  trial  value  of  t  (instead  of  y/z)  must  now  he  found 
which  satisfies  equation  (Cl) . 


5T 
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TABLE  I 

SEPARATION  POINTS  CALCULATED  BY  APPROXIMATE  METHOD  (EQUATIONS  (25)  AND  (26)) 
COMPARED  WITH  STEWARTS  ON’S  SOLUTION  AND  WITH  OTHER  CALCULATIONS 


TABLE  II 


VALUES  OF  CRITICAL  REYNOLDS  NUMBER  FOR  STABILITY  AS  CALCULATED 
FROM  APPROXIMATE  SOLUTION  (EQUATION  (22))  OF  EQUATION  (l4) 
COMPARED  WITH  OTHER  SOLUTIONS  FOR  ZERO  MACH  NUMBER 
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TABLE  III 

COMPARISON  OF  CALCULATED  SKIN-FRICTION  COEFFICIENTS 


Method 

Roo^/^Cf/l 

Fourth-degree  profile  (equation  (35)) 

2.390 

Sixth-degree  profile  (equation  (32)) 

2.530 

Solution  of  Falkner  and  Skan  (reference  l4) 

2.468 

TABLE  IV 

COMPARISON  OF  CALCULATED  CRITICAL  REYNOLDS  NUMBERS  FOR  STABILITY 


Method 

Foirrth-degree  profile  (equation  (35)) 

Sixth-degree  profile  (equation  (32)) 

Solution  of  Falkner  and  Skan  (reference  l4) 

210  X  106 

170 

240 

4o 
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Figure  4.-  Comparison  of  fourth-  and  sixth-degree  velocity  profiles 

with  solution  of  Falkner  and  Skan  (reference  l4).  uj^/uoo  =  I;  =  0, 
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0  .4  .8  1.2  1.6  2.0 


X  /iiooL 

L  V  CV^ 

Figure  7*-  Velocity  profiles  based  on  approximate  solution  (equation  (36)) 
of  ordinary  differential  equation  (equation  (l4));  foiirth-degree 
velocity  profile,  ^±/^cx>  =  I  =  OA. 
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